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SUMMARY 


The accuracy of the method, presented, in NACA Rep. No. 824 for rapidly 
calculating the increment in the velocity distribution about a symmetrical 
airfoil due to angle of attack has been investigated by use of the exact 
potential-f low-theory solution of the velocity distribution about an 
arbitrary airfoil given in NACA Rep. No. 452. The analysis shows that the 
rapid method described in NACA Rep. No. 824 evaluates the effect of angle 
of attack on the net velocity distribution at any point on a symmetrical 
airfoil with a maximum error which can be expressed as an increment in 
velocity ratio by the relation (l - cos a) where a is the angle 

of attack and \fS^ is the velocity ratio on the airfoil at zero lift 
calculated according to the exact theory. 


INTRODUCTION 


A method Is described in reference 1 which permits the rapid esti- 
mation of the pressure distribution at any lift coefficient about large 
numbers of airfoils from a rather limited amount of theoretical data. 

The method is. based on the assumption that the velocity at any point on 
an airfoil can be resolved into components which are attributable to the 
induced velocity on the basic thickness form at zero lift, the load 
distribution due to angle of attack, and the load distribution due to 
camber. The accuracy of the method described in reference 1 for obtaining 
and applying the increment in velocity associated with the loading due 
to angle of attack has been open to some question. The purpose of the 
present paper is to give an exact account of the method by means of which 
the Increment in velocity due to angle of attack is obtained and the 
inaccuracies to which the assumptions involved in the method may lead. 

The analysis is concerned only with symmetrical airfoils and does not 
consider any possible effects of camber on the velocity increments due 
to angle of attack. 


2 


NACA TN 1884 


COEFFICIENTS AND SYMBOLS 


X 


n 

a 

S 

Ho 

P 

*0 

\/S 


V 

At. 


t'R 


4 

K 

*oi 


reference length 
airfoil chord 

distance along airfoil chord 
section lift coefficient 
section normal-force coefficient 
section angle of attack 

pressure coefficient 


(hr) 


free-stream total pressure 
local static pressure 
free-stream dynamic pressure 
velocity ratio for incampressihle flow 
local velocity at some point on the airfoil 
free-stream velocity 

increment in velocity ratio caused ty additional type loading 
due to angle of attack 

resultant, pressure coefficient} difference between local 
upper-surface and lower-surface pressure coefficients at 
the same chordwise position 

circulation 

parameters which are functions of airfoil geometry (reference 2) 


Subscripts; 

U upper surface 

L lower surface 

o conditions calculated for airfoil at zero lift 
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DISCUSSION 


The equation given in reference 1 for obtaining the net velocity 
ratio v/s at some point on a symmetrical airfoil at a lift coeffi- 
cient ci is 

(i) 

Values of \/s^, the velocity ratio for the symmetrical airfoil at zero 
Av a 

lift, and 1 ^ , the incremental velocity ratio due to angle of attack for 

a lift coefficient of 1-0, are tabulated in reference 1 for a large 
number of airfoils- The sign convention followed is that all velocities 
which cause a flow in the direction of the trailing edge are positive. 

AVq 

A negative sign is therefore employed with —— when the net velocity 

ratio VS is calculated for the lower surface at positive lift coeffi- 
cients. 

The values of y/s^ presented in reference 1 for the various basic 
thickness forms at zero lift were calculated by the exact theory of 

Av a 

Theodorsen and Garrick (reference 2) . The values of were calcu- 

lated by means of the following relation*. 


Av a = VS " JEL (2) 

V 2 

where \fS\j and ^/s £ are the upper-surface and lower-surface velocity 

ratios, respectively, calculated by the exact theory for some lift coef— 

Av a 

ficient c^ which is usually small. Values of for a lift coef- 

ficient of 1.0 were obtained by linearly scaling the values given by 
equation (2). Two basic assumptions are implicit in equations (l) 
and (2): first, that the mean of the upper- and lower-surface velocity 

ratios at any lift coefficient, given by 

fzz - & ; & < 3 ) 


is the same as the velocity ratio calculated for the same point on the 
airfoil at zero lift and second, that the incremental velocity 
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ratio is a linear function of lift coefficient. The validity of 

these two assumptions will now he investigated. 

From the exact potential-flow theory (reference 2), the following 
equations may he written for the velocity ratios on the upper and lower 
surfaces of a symmetrical airfoil: 


Jsjj = K[sin(a + 0) + sin a| 

(4) 

= K[sin(a. - 0) + sin a] 

(5) 


where K and 0 are functions of only the airfoil geometry and have the 
same value at corresponding points on the upper and lower surfaces of a 
symmetrical airfoil. The angle of attack is given by a. The velocity 
distribution about the airfoil at zero lift is, of course, given by the 
equation: 

^ = K Bin 0 ( 6) 


Since 0 is negative on the lower surface, equations (5) and (6) show 
that velocities which cause the flow to progress toward the trailing 
edge on the lower surface will have a negative sign. In order to be 
consistent with the previously defined sign convention, however, a nega- 
tive sign must be placed in front of equation (5). Thus 

y/s£ = “K[sin(a - 0) + sin a] (7) 


An indication of the accuracy of the assumption that 


\/So - \/Smean 


caii be obtained by substituting the expressions given by equations (4) 
and (7) in equation ( 3 ) and determining the manner in which the resultant 
expression differs from the exact equation (6) for >/s^. Performing 
the suggested substitution 


n/®U + •_ r 

r " 

[sin(a + 0) + sin aj - [sin(a - 0) + sin a] 

2 . 

2 ” 
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and expanding and simplifying give the following expression 


^ ^ = VSmean = K-sin 0 cos a 


( 8 ) 


Equation (8) for the mean velocity ratio \/Smean is seen to differ from 

the exact expression for the velocity ratio about the symmetrical airfoil 
at zero lift (given by equation (6)) only by the factor cos a. Since 
cos a is nearly 1.0 for relatively small values of a, such as would be 

of interest, the assumption that = v/Smean is seen to be justi- 

fiable. Furthermore, the approximation is seen to be in no way dependent 
upon the chordwise position along the airfoil. 


In order to investigate the assumption that 



varies linearly 


with lift coefficient, the expressions given by equations (4) and (7) are 
substituted in equation (2) for the incremental velocity ratio 


Av a 

S — ir^ 

[sin(a + 0) + sin aj + 

sin(a - 0) + sin d] 

_ V 

v " K 

2 




(9) 


where the upper-surface and lower-surface velocity ratios given by 
equations (4) and (7) are for some lift coefficient Simplifying 

equation (9) and dividing both sides by the lift coefficient c^ yield 


^a = K sin a / CQa 0 + £\ (10 ) 

Yc z c z V J 

An expression relating sin a and cj can be obtained from the following 
equation for the circulation F taken from reference 2 

\|f 

T = 4rtVae °sin a (ll) 


where is a constant for any given airfoil, V is the free-stream 

velocity, and a is a constant which is approximately one-quarter of 
the airfoil chord c. The relation between a and c varies slightly 
for different airfoils, but this variation has no bearing on the present 
discussion. The lift coefficient and the circulation are related by the 
following expression: 
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Substituting from equation (12) in equation (ll) gives the following 
relation between sin a and cj 

sin a = (13) 
2ite 

/ s . . A Vo, 

The use of equation (13) with equation (10) for gives 


Ava 

V 


Kjcos 0 + l] C 2 


2ne 


(14) 


Equation (l4) shows that the incremental velocity ratio does, in fact, 
vary in an exactly linear manner with lift coefficient. Equation (14) 
also provides a very convenient means for calculating the values 

of since all the parameters involved in equation (14) must be 

evaluated for the calculation of the velocity distribution about the 
airfoil at zero lift. 


The analysis presented has shown that for symmetrical airfoils the 
mean of the exactly calculated upper-surface and lower-surface velocity 
ratios for any lift coefficient is related to the exactly calculated 
velocity ratio at the same point on the airfoil at zero lift by the 
following expression: 

x/Smean = Vs^ cos a (15) 


and that the incremental velocity ratio varies in an exactly linear 

manner with lift coefficient. The only error which will be incurred in 
the net velocity ratio as a result of employing equation (l) together with 
the values of calculated for the airfoil at zero lift and the values 

Av 

of calculated from equation (2) or (14) is therefore an increment in 

velocity ratio on the upper and lower surfaces which is given by (l - cos ct)\js ~. 


With the aid of the results just obtained, the relation is now shown 
between the lift coefficient c\ employed in equation (l) , for which the 
pressure or velocity distribution is desired, and the lift coefficient which 
is obtained from an integration of the pressure distribution calculated 
according to equation (l) . The resultant pressure coefficient at some 
point along the airfoil chord Pp> is defined as the difference between 
the upper-surface and lower-surface pressure coefficients at that point. 

An integration of P E along the chord of an airfoil at some angle of 
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attack gives an expression for the normal-force coefficient c n 


P1.0 



(16) 


The lift coefficient is then given hy the relation 


c l 


cos a 


ni.o 


r R 



(17) 


since the drag is zero for perfect fluid, flow. In terms of the exactly 
calculated values of the velocity ratio for the basic thickness form at 

- Av«a 

zero lift \JS Q and. the values of for the lift coefficient under 

consideration, the resultant pressure coefficient is given approximately 
by equation (l) as 



The substitution of the approximate expression for Pp (equation (l8)) in 
equation (l6) gives 


c n ~ 


P1.0 


0 


V \o, 


(19) 


whereas the accurate expression for c n would be 


Pl.O 


c n _ 


^ \/^me 


an 


Av a p /x 


V Vc 


( 20 ) 
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_ \J^m> 


Since Js^ = equation (19) can be written 

C OS CL 


ni.o 


p ^ — — — 

11 cos a 


UO 


,VS= V s <f) 


Thus, the use of the approximation v/s^ « v/Smean is found, to result in 
the evaluation of the lift coefficient cj instead of the normal-force 

p 

coefficient c n , inasmuch as C 7 = . 

11 1 cos a 


CONCLUDING REMARKS 


The accuracy of the method described in reference 1 for rapidly 
calculating the increments in the velocity distribution about a symmet- 
rical airfoil due to angle of attack is in error only by a factor cos a 
over the entire chord- In view of the discrepancies, resulting from 
the effects of viscosity, which are known to exist between airfoil 
velocity distributions calculated by any perfect fluid theory and those 
obtained experimentally, it is difficult to see how any situation might 
arise in which a theoretical, perfect-fluid velocity distribution about 
a symmetrical airfoil would be required to a degree of accuracy greater 
than that provided by the method described in NACA Rep. No. 824. 
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